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Abstract 

In recent works [jlj] (BHP), a generalized universality has been proposed, 
linking phenomena as dissimilar as 2D magnetism and turbulence. To test 
these ideas, we performed a MC study of the 2DXY-m.ode\. We found that 
the shape of the probability distribution function for the magnetization M is 
non Gaussian and independent of the system size -in the range of the lattice 
sizes studied- below the Kosterlitz-Thoules temperature. However, the shape 
of these distributions does depend on the temperature, contrarily to the BHP's 
claim. This behavior is successfully explained by using an extended finite-size 
scaling analysis and the existence of bounds for M. 
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The study of critical phenomena is of great interest not only because it allows the un- 
derstanding of a large number of very different physical systems, like the super fluid Helium 
three, low temperature superconductors, ferromagnetic-paramagnetic systems, turbulent flu- 
ids and plasmas, polymers, snow flakes and earthquakes, but also due to the existence of 
scale independence of the fluctuations at the critical temperature. In fact, although the 
underlying inter-molecular forces, responsible for the existence of phase transitions, have a 
well-defined length scale, the structures they give rise do not. This leads, very close to the 
critical temperature, to the power-law behavior of physical quantities, which characterizes 
universality. The main challenge of the theory of critical phenomena is to explain how dis- 
similar systems exhibit the same critical behavior. Renewed interest in this subject has been 
raised, because in a seminal paper || (BHP), it was argued that turbulence experiments 
can be explained in terms of a self-similar structure of fluctuations, just as in a finite crit- 
ical system like the harmonic finite 2DXy-model {2D HXY -model). The starting point of 
this conjecture was the observation that the probability distribution function (PDF) of the 
injected power fluctuations in a confined shear turbulent flow |2| has the same shape as the 
PDF of the magnetization in the 2DHXY-model. It was also proposed that this analogy 
should provide a new application of finite size scaling in critical systems with experimental 
consequences. 

In this paper, we report the results of a high precision Monte Carlo study of the full 
2DXY-model. This computation was carried out over the whole physical range of tempera- 
tures. The magnetic susceptibility was computed and the lattice-shifted critical temperature 
was obtained for different lattice sizes. Scaling laws for the magnetization-temperature ratio 
were tested. Our results agree with the rigorous findings of Chung ||. We also found that, 
below the Berezinskii-Kosterlitz-Thouless temperature Tbkt, the shape of the PDF of the 
magnetization is non-Gaussian and independent of the lattice size, in agreement with pre- 
vious results ||]. However, we found that the shape of these distributions do depend on the 
temperature, contrarily to the generalized universality claimed in (BHP), who stated that 
the PDF of magnetization is independent of both, system size and temperature. Our results 
can be seen as a powerful extension of finite-size scaling and phenomenological renormaliza- 
tion of the PDF, suggested originally by Binder in the context of the Ising model J5|, with 
sligth modifications introduced by the bounds of M. This allows in particular to understand 
the scaling form of the PDF of the order parameter in the 2DXY-mode\ as well as in the 
turbulent system. 

The universality proposed by (BHP) might go beyond the idea of equivalence classes 
in Wilson's renormalization group approach ||, by including into a generalized universality 
class systems sharing the properties of finite size, strong correlations and self-similarity, even 
if their space-dimensions are different. 

In the two-dimensional probability distribution for the magnetization is calculated 
by means of a Monte Carlo simulation in the context of the 2DHXY-model. This model is 
a further simplification of the Villain model f7j, where the vortex variable n is not a thermo- 
dynamical quantity, but it is constrained to the values n = —1, 0, 1. By using diagrammatic 
techniques, they showed that this asymmetry could be the result of three-spin interactions 
and higher order corrections. 

Here, we consider the 2DXY-model, which describes classical planar spins with nearest 
neighbor interactions, with a Hamiltonian given by 



2 



H=-JY,cos(6 i -e j ) (1) 

where J is the ferromagnetic coupling constant and d\ is the angle of orientation of the unitary 
spin vector ifj. The summation is over nearest neighbors and the spins are defined 
on the sites of a square lattice of lattice size L, with periodic boundary conditions. From 
hereon the ratio k B /J is set equal to unity throughout the paper. This model undergoes a 
remarkable binding-unbinding topological phase transition, such that the free energy and all 
its derivatives remain continuous ||, and no long-range order at low temperatures exists, as 
stated by the Mermin- Wagner theorem || . This model has been extensively studied through 
both numerical and analytical methods |TD . 



Our simulation was performed on a square lattice of lattice sizes L=10, 12, 16, 22, and 
32 respectively. We estimate the MC sweeps needed for thermalization by plotting some 
observables like magnetization and energy. Typically 10 5 MC sweeps were used to reach 
thermal equilibrium. For thermal averages we used 5 x 10 5 spin configurations otj. Because 
the 2DXY-mode\ has a continuous line of critical points below Tbkt-, special care was taken 
to choose statistical independent configurations to evaluate thermal averages of physical 
observables X. This was achieved by computing its normalized autocorrelation function 
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nl ts\ < X ai X ai+K > - < X ai >< X aK > 

C{K) = <xi>-<x a , 



(2) 



where X a . is the value of X in the configuration a« at the i — th step along the MC-path 
through the configuration space, and the average < . . . > was taken over this particular 
path of configurations separated by K steps from each other. C{K) — 1 for K = 0, but for 
large enough K, C(K) drops to zero, which means that these configurations become totally 
uncorrelated. We choose K so that C(K) was less that the recommended value 0.05 |12| . It 



is well known that as a critical system approaches the critical temperature, the decorrelation 
time t diverges with the power law r ~ £ 2 , where £ is the (divergent) correlation length of 
the system and z is known as the dynamical critical exponent, which is approximately two 
for local-flip algorithms like the Metropolis algorithm. This phenomenon is known as critical 



slowing down |13fl . This means in practice a serious limitation to numerical simulations of 
critical systems close to a critical point. 

In fig. 1 we show MC data for the susceptibility for I? = 256 spins, as a function of 
the temperature. The peak occurs at the value Tc(16) = 1.15, and corresponds to the 
temperature at which the correlation length equals L, which is the standard definition of the 
critical temperature of a finite system. We compute also the errors (standard deviations), 
which become larger as the critical temperature is approached. Another interesting feature 
of these errors is that they are larger below Tc- This can probably be explained because of 
the comparative larger correlations lengths in this region, which corresponds to a continuous 
line of critical points with temperature dependent exponents in the infinite volume limit. 

We computed the critical temperature for the lattice sizes L=10, 12, 16, 22, 32. For 
L=32 we found an effective transition temperature Tc = 1.08, in agreement with the value 
obtained by [[TJJ, where the linearized RG equations for the finite size scaling were used. 
Fig. 2 shows Tc(L) as a function of (ln(L)) -2 . The values can be described by the finite 
size scaling formula || 
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where is the extrapolated value of the critical temperature for infinite volume. Within 
a few percents of error, we found that the value of agrees with the seemingly exact 
value 0.892 of the critical temperature Tbkt of the Berezinskii-Kosterlitz-Thouless phase 
transition. 

The ratio of the mean magnetization to critical temperature is plotted in fig. 3 as a 
function of ln(L). These values are compatible with a negative straight line, as suggested 
by 0] in the context of the harmonic XY-model. The values closer to the origin have larger 
statistical errors probably because of the finite size effects, which are proportional to ln(L) 

IJH- 

Finally it should be emphasized that we did not make use of scaling relations to define 
physical quantities, like the shifted BKT-temperature T*(L), or temperature at which the 
renormalized spin- wave stiffness becomes the universal value 2/ir of the infinite system. 
This is because the use of the BKT theory beyond its confirmed validity needs at least 
justification, (the scaling region is defined by the inequality | T — Tbkt) /Tbkt |< 1CT 2 , 
where the renormalization group equations confidently apply fl7l). In fact, we were able 
to obtain accurate values for thermal averages and test the scaling equations (see eqn. @ 
and Figs. 2 and 3), in spite of the difficulties of numerical simulations due to the essential 
finite-size effects present in this model || (ln(L) dependence of physical quantities), and the 
very narrow critical region. 

The renewed interest in the PDF of fluctuations of magnetization, M, is a consequence 
of the observation that similarly shaped PDF arise in completely different systems. For 
instance, in || it was found that fluctuations of the injected power in a confined turbulent 
flow show the same behavior. In Fig. 4 a plot of aQ(M) as a function of (M— < M >)/er 
can be seen, for lattice sizes L = 16 and L = 32 at the same temperature T = 0.70. Here, Q 
is the PDF of M and corresponds to Pl{M) in the language of Binder (see discussion below 
eqn. (|J)). These curves have similar shapes like those found in turbulence experiments, but 
only within a reduced range of temperatures below Tbkt- These PDFs can be conveniently 
compared with the universal form Tl(y ) = K(e b( ' y ^ s ^ eb( ' v s) ) proposed in ||T|, by plotting the 
ratio crQ(M; T)/H(M) vs. (M — (M))/a. In Fig. 5, four of such plots are displayed for 
T = 0.40,0.80,0.95 and 1.05. The upper curves are successively multiplied by factors of 10 
for clarity. As can be seen, when T is increased, these ratios consistently change, showing 
the dependence of Q(M; T) on the temperature. The rightmost part of the lower curves is 
raised, probably due in part to the upper bound M — 1. On the other hand, for T j = 0.95 we 
can see that the opposite occurs. This temperature is slightly greater than T*, above which 
the population of spin vortex pairs begins to increase. This happens because the system 
volume occupied by these vortices no longer contributes to the magnetization, which leads to 
a depleted probability density. This can be appreciated in the central and rightmost parts of 
the two upper curves. The hills are due to the bounds in the magnetization < M < 1. At 
T = 1.05 this effect is greatly enhanced, and the leftmost part of the curve shows even more 
clearly the effect of the lower bound M = 0. Concerning turbulent flows, we do not expect 
this type of bounding effects in the statistics of injected power. In principle, there are no 
limits to the fluctuations of such a quantity, and negative values are not excluded, meaning 
that the flow is delivering power to the driving system. Although this type of events are 
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expected to be very unlike, they are not forbidden. 

The temperature dependence of the PDF is not a surprising result. In fact, the use of the 
probability distribution of the order parameter to study finite size scaling and phenomeno- 
logical renormalization, has been discussed by Binder in the context of the Ising model 



For the region £ ~ L he proposed that the probability distribution function Pl(M) does not 
depend separately on the three variables £, L, M, but only on the two scaled combinations, 
L/t and MC ": 

P L (M) = ^ lu P(L /£, M^ /v ) = L p/v P{L /£, ML p/v ). (4) 

He also argued that in the critical region £ 3> L, Pl(M) is no longer Gaussian. In the 
scaling region, it is a good approximation to take P L (M) equal to the PDF proposed by 
Bramwell et al. (The standard deviation a plays the role of L~^l v in the BHP distribution; 
this can be seen by using the standard definitions of the critical exponents and the relation 



a — y jpX )• Nevertheless, and away from the region defined by | T — Tbkt \ /Tbkt < 1CT 2 , 
there is a temperature dependence in expresion (4) via the correlation length £ for finite 
size. It turned out that this dependence is rather weak in the range 0.5 < T < T*(L) 
(T*(L = 16) ~ 0.94), but out of this range this dependence becomes stronger due to the 
presence of vortices and/or bounds. 

In conclusion, we found that the probability distribution function for the magnetization 
is indeed independent of the system size, but its shape happens to vary with the temperature 
of the system, contrarily to the generalized universality proposed by Bramwell et al. . This 
effect comes from the intrinsic temperature dependence on the first scaled variable L/£ in 
the distribution function proposed by Binder for the order parameter M. Also, there is a 
contribution coming from the constrained character of the magnetization. 

We thank L. Vergara for helpful suggestions. This work was supported in part through 
projects FONDECYT 1980608 and 1990169, and DICYT 04-9631PA and 04-9631LM. 
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Figure captions 

Figure l:Susceptibility for L = 16 in the range < T < 3. The peak at 1.15 corresponds 
to the shifted critical temperature. 

Figure 2: The shifted critical temperature for different lattice sizes is plotted as a 
function of the system size. 

Figure 3: Scaling relation for the magnetization-temperature ratio as a function of 
system size. 

Figure 4: Plots of aQ(M) vs. (M— < M >)/a at T = 0.70 for lattice sizes L = 16 (+) 
and L = 32 (*). 

Figure 5: aQ(M) /TI(M) ratios for four values of temperature (see text). 
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